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Abstract 

We introduce the idea of natural mass matrices, an organizing 
principle useful in the search for GUT scale quark mass matrix pat- 
terns that are consistent with known CKM constraints and quark mass 
eigenvalues. An application of this idea is made in the context of SUSY 
GUTs and some potentially "successful" GUT scale mass patterns are 
found. The CKM predictions of these patterns are presented and some 
relevant strong CP issues are discussed. 
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1 Introduction 

It has been a theoretical quest for nearly 20 years to devise interesting mass 
matrix patterns which could provide sound predictions for the CKM ma- 
trix and quark mass eigenvalues Recently, these efforts have centered 
on constructing patterns at the GUT scale in SUSY theories J3|. Most 
of these attempts, although quite successful, have failed to produce results 
in complete agreement with precise low energy data(LED)Q. Part of the 
difficulty lies in the fact that one has to rely upon "guesses" at the GUT 
scale which are then extrapolated down to the weak scale, with the hope 
that the mass matrices so obtained will give rise to acceptable fits to the 
LED. A somewhat more promising approach perhaps, would be to reverse 
the process: constructing LED-consistent mass matrices at some low energy 
scale and then evolve them upwards to see if interesting GUT patterns would 
emerge. However, since there exists considerable arbitrariness in this con- 
struction, one typically has to be content with studing only certain special 
cases H || often guided merely by simplicity. Here, we would like to suggest 
an organizing principle which may be helpful. This is the idea of natural 
mass matrices, which severely restricts the aforementioned arbitrariness in 
the mass pattern construction, thereby allowing a search for viable GUT 
patterns more systematically and efficiently. 

This paper is organized as follows: in Sec. 2, we summarize our present 
knowledge of the CKM parameters and quark masses; in Sec. 3, we introduce 
the idea of natural mass matrices along with our method of mass-matrix- 
parametrization that facilitates its implementation; in Sec. 4, we present 
some Hermitian GUT scale mass patterns that emerge from this approach 
and their CKM predictions; in Sec. 5, we discuss issues connected with Her- 
miticity breakdown and associated problems with strong CP violation; finally 
in Sec. 6, we present our conclusions. 

2 CKM and Quark Masses: a Summary of Present Status 

2.1. CKM Matrix 

In its standard form [jTJ, the Cabibbo-Kakayashi-Maskawa (CKM) matrix 
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IS 

cic 3 sic 3 s 3 e~ lS 

[CKM] S = I -S1C2 - ciS2S 3 e 115 dc 2 - sis 2 s 3 e i<5 s 2 c 3 j . (2.1) 

S1S2 - cic 2 s 3 e J<5 -C1S2 - siC2S 3 e 4<5 c 2 c 3 

Taking into account of the experimental hierarchy in the mixing angles, one 
can write 



si = sin#! = A ~ 0.22 , 
s 2 = sin 82 = AA 2 , 

s 3 = sin ^ 3 = AaX 3 , (2.2) 

with A, a being parameters roughly of 0(1). The CKM matrix then takes 
the Wolfenstein |Sj form Q 

/ 1 _ 1 A 2 - |A 4 A Aa\ 3 e-* 5 \ 

[CKM] S = -A 1 - ±A 2 - (±,4 2 + |)A 4 ,4A 2 

\ AX 3 (l - ae iS ) -AX 2 + ±AA 4 1 - ±A 2 A 4 J 



+ 0(A 5 ). (2.3) 

The parameter A is fixed by V cb which, in turn, is determined from 
semileptonic B decays. The most recent analysis from CLEO data || gives, 

I = 0.0378 ±0.0026 A = 0.78 ± 0.05 . (2.4) 

The parameter a is fixed by the ratio |\4b|/|Vd>|. This in turn can be 
extracted from a study of semileptonic B decays near the end point region 
of the lepton spectrum, where b —>■ c quark transitions are forbidden. The 
most recent analysis H gives, 

\V ub \/\V cb \ = 0.08 ±0.02 a = 0.36 ±0.09 . (2.5) 

The phase 5 (or the parameter 77) can be gotten by combining the mea- 
surements of the e parameter in K — K mixing and those of ArriB d in — 

1 Wolfenstein uses the parameters p, r/ instead of a, 8. They are related by ae~ lS = p—irj. 
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mixing with the value of |14&|/|Vd,|. A recent analysis |T(J gives an allowed 
region in the p — r\ plane roughly specified by the ranges 

r] ~ [0.2, 0.5] , p ~ [-0.4, 0.3] (2.6) 

with a corresponding CKM phase 

5 ~ [45°, 158°] . (2.7) 

Quark Masses 



For the purpose of calculating quark mass ratios and their RG scaling, it 
is convenient to express all quark masses as running masses at some common 
energy scale. We shall choose this scale here to be the mass of the top 
quark m t . We summarize below what is known about the quark masses and 
then extrapolate all the results to the scale m t . From the recent "discovery" 
papers on the top quark fill , one infers a value for the physical mass mf hys , 
which is related to the running mass by 



m t {m t ) 



ni 



Phys 



1 + ^a s (m t ) 



These results [11] suggest that 



m t (m t ) ~ (175 ± 15) GeV. 

For medium heavy quarks, the analyses of charmonium and bottomonium 
spectra [12| give 



m c (m c ) = (1.27 ±0.05) GeV , 
m b (m b ) = (4.25 ±0.10) GeV . 



Finally, for light quarks, current algebra analyses |12] give the following values 
for the masses at a scale of lGeV, 



m u {\GeV) 
m d (lGeV) 
m s (lGeV) 



(5.1 ± 1.5) MeV 
(8.9 ± 2.6) MeV 
(175 ± 55) MeV 
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These mass values are individually uncertain to 0(30%) but are rather better 
constrained [13] [14] by the current algebra relation 



(— ) 2 + i (— ) 2 = 1 . with Q = 24 ± 1.6 • (2.8) 
m d Q 2 m d 

The RG scaling of the medium heavy and light quark masses to m t are 
calculated to 3-loops in Ref. |§, with the result sensitive to the precise value 
of the strong coupling constant a s (M z ). Using a s (M z ) = 0.117 ± 0.05 and 
expressing all quark mass ratios in terms of the small parameter A ~ 0.22, 
the above results allow one to write the diagonal quark mass matrices as 



m t 



175 GeV 



( iuA 1 o o 

£ rf A 4 
V i tt 



= m t {m t ) Mf" 9 , 
where £ ut = 0.49 ± 0.15, U = 1.46 ± 0.13, £ u = 1 ± 0.09; and 

UA 4 

U^ 2 

£ bb 
= m b (m t ) Mf a9 , 
where ^ = 0.58 ± 0.18, ^ = 0.55 ±0.18, £ bb = 1 ± 0.05. 



Mf ag {m t ) = 2.78 GeV I £ sb \ 2 



(2.9) 



(2.10) 



3 Natural Mass Matrices 

3.1. A Heuristic Two-Generation Example - the Notion of Naturalness 

To introduce the idea of natural mass matrices, we consider first the 
simple two-quark-generation case, assuming that these mass matrices are 
Hermitian. A general 2x2 Hermitian mass matrix for the first two quark 
families can always be rewritten, after some trivial phase redefinitions of the 
quark fields, as some real symmetric matrix. Two such matrices for the u- 
quarks and d-quarks: M u = m c M u , Md = m s Md, in turn can be diagonalized 
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by some orthogonal matrices O u and Od, resulting in a Cabbibo-quark-mixing 
matrix C. One hasf] 

T u M u O u = M^ = 

T d M d O d = M d d iag = ( sas n ' ; ) (3.1) 
and 

C = OlO d ={ ™ B ° M A . (3.2) 
^ - sin 9 C cos 9 C J 

In the above, A = sin Q c as before and the matrices O u , Od have the same 
form as the matrix C with angles 9 U , 9 d instead of 9c- We notice that it 
follows from Eq.( |3.2| ) that 9c = Od — 9 U . Moreover, since the matrix C is 
invariant under the changes 

o d ^oo d - o u ^oo u , 

where O is some arbitrary orthogonal matrix, we see that M u and Md are 
fixed only up to a common similarity transformation 

M u «-> T M U O; M d ^ T M d O - 

Because of this freedom and since 9c ~ A < < 1 , we can always arrange to 
have both 9d « 1 and 9 U « 1. PI We can now contemplate three different 
options for the angles 9 U and 9d- 

(i) sin 9d ~ A , sin 9 U ~ A ; 

(ii) sin0 d ~A, sin6> u <A 2 ; (3.3) 
(Hi) sin 6*^ < A 2 , sin 6* M ~ A . 

Upon examining the expressions for the mass matrices M u and Md- 



M u u ^ u 



n -Cjrdiagr^T ^ ( £«cA + sin 9 U sin 9 U \ 

Uu u Uu ~ { sin9 u 1 J ' 

OdMf^Ol * ( UX ] + ff ° d Sin / d ) , (3.4) 



2 Note that £„ c = WfctA) = (1.53 ± 0.49) while fa = fa/fa = (1.05 ± 0.47). 

3 As will become apparent below, this restriction need not to be imposed separately 
for it is in conformity with our naturalness requirement. It is done here for ease in the 
discussion that follows. 
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one sees that if options (i) and (Hi) were to hold, one would require a severe 
fine-tuning of the matrix element [MJ11, forcing [M u ]n ~ ([Af u ]i 2 ) 2 to arrive 
at the large m u /m c ~ A 4 experimental hierarchy. Such a fine-tuning appears 
to be unnatural. For the two-generation case, the natural mass matrices have 
the form 

M„ , ( °f );*=( f ) , (3,) 

corresponding to option (ii) in which 



UC 1 



sin 9 U = a u \ 2 , a' u - ct u = £ 
sin 9 d = a d \ , a' d - aj = 

and, 

a d - a u X ~ 1 . 

With the parameters a and a' of 0(1), one gets the observed hierarchy (i.e. 
the £'s being of 0(1)) without any need for fine-tunings. 



3.2. Three- Generation Extension 



It is straightforward, though somewhat tedious, to extend the idea of 
natural mass matrices to the three-generation case. To facilitate its imple- 
mentation however, we need to introduce a convenient parametrization of the 
mass matrices based on a perturbative expansion in "A" . The procedure we 
shall adopt is analogous to, but a slight generalization of, a method employed 
by Ramond, Roberts and Ross 0. The benefits of our generalization, aside 
from allowing a search for natural mass patterns, also include the flexibility 
of simultaneous adjustments of matrix elements in M u and M d (useful in 
imposing "0"s or arranging for equalities among them). 

Consider some general 3x3 Hermitian mass matrices M u = m t (m t )M u 
and M d = mb(m t )M d . These matrices can be diagonalized by some unitary 
matrices U and D, resulting in a CKM-quark-mixing matrix [CKM]: 

M u = UMt ia9 U\ (3.6) 
M d = DM d iag D\ (3.7) 
[CKM] = U ] D . (3.8) 
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To proceed, it is useful to make two observations: 

(1) As in the two-generation case, if we change U — > NU and D — > ND 
(where N is some arbitrary unitary matrix), the matrix [CKM] remains 
unchanged. The matrices M u , M d are thus unique up to an arbitrary (but 
common) unitary transformation: M u <-> N'(M U )N, Md <-> N'(M d )N. As a 
result of [CifM] ~ 1, this arbitrariness allows us to focus only on "small" 
transformations, i.e. U ~ 1, D ~ 1. Furthermore, in the particular case 
where iV = 0^, 0^ being some phase matrix, the induced changes in M u and 
Md {M u — > <pLM u (p^ L ; Ma — > (jy^Md^i) can be absorbed by a trivial quark 
field phase redefinition. 

(2) So far as constructions of the mass matrices M u and Md out of the 
matrices Mf ag , Mf ag and [CKM] are concerned, a change of [CKM] -> 
0t [CiTM] 0d <^=^> D — > Dcf) d ; ^ — > Ucf) u , with the 0's some arbitrary phase 
matrices, is inconsequential due to the fact that (p u M^ a9 0^ = M^ ag and 
<j) d Mf ag <j)l = Mf ag . 

Based on these observations, one is always justified working with a specific 
form of the CKM matrix. Thus one can write 

U^D — [CKM] = <pl [CKM] s <f)d 

with 

D = 4>{D s (t> d ■ U = ^ L U s <j) u . 

The phase matrices 0^, <pd can be further chosen so as to render D s to 
have the same form as the matrix [CKM] S of Eq. (|2.1| ), and it follows that 
U\D S = [CKM] S . (Notice that having chosen D s to be of the same form as 
the matrix [CKM] S , there is now no more freedom to redefine U s and U s , in 
fact, will not have quite the standard CKM form.) The mass matrices M u , 
M d constructed as 

M u = U s Mt ia9 U\- 

M d = D s Mf a9 D\ (3.9) 

are still perfectly general for our considerations.^ Having established their 
generality, we shall drop the subscript "s" in the matrices [CKM] S , D s and 
U s hereafter for notational brevity. 

4 The matrices M u and Md in Eq.( |3.9| ) can still be changed by an overall common 
unitary transformation. It is possible to remove this remaining freedom by adopting a 
(slightly) unconventional parametrization of the CKM matrix. Defining M u and Md as 
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With these preliminaries out of the way, it is straightforward to give 
explicit parametrizations for the matrices U and D. For these purposes, it is 
useful to note that, in its standard form, the matrix [CKM] can be expressed 
as 

[CKM] = C 2 AC 3 A^d 




with 



A 




Because D, by assumption, also takes the CKM form, we can likewise write 

D = C 2d A d C 3d A ] d C ld . (3.10) 

Here the matrices C id and A d are defined analogously to C and A, except 
they involve some new angles 6 id (i = 1,2,3) and 5 d . It is convenient, in 
addition, to define three more orthogonal matrices C{ U (of the same form as 
the Cj's) involving angles 6 iu {i = 1,2,3) and satisfying the relations 

CjvPid = C% (no sum over the index "i" ) . 



In analogy to the two-generation case, the angles 9 iu , 9 id and 9i obey 



0i — 9 id — 9i 



1,2,3) 



(3.11) 



Using these definitions, the matrix U is easily seen to be 

U = D[CKM} ] 

= {C 2u }{C 2 (A d C 3u A d )Cl} 

{C 2 (A d C 3 A d )C lu (A Cj At) C|} 



(3.12) 



in Eqs.(p.6[) and (3/7), one has always the freedom to choose D to be of the standard 
CKM form D = D s , so that the matrix Ma is described by 6 real variables and 1 phase. 
However, then there is no further freedom of redefinition, hence U being a general 3x3 
Hermitian matrix must involve 6 real variables and 3 phases. It is easy to show that 
this more general parametrization corresponds to a definition of the CKM matrix (in the 
standard form) [CKM] S = U^D S $, with $ being a diagonal phase matrix containing 2 
arbitrary phases. We find it more convenient to take $ = 1 and have the matrices M u 
and Md undetermined by an overall common unitary transformation. 
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This matrix is not quite of the standard CKM form. However, since the ma- 
trices Cj's ~ 1, the matrix U is not that different, apart from the placement 
of some phase factors. 

Once one has explicit forms of the matrices U and D, the "A" expansion 
of these matrices and hence of M u , Md is accomplished by letting: 

e ld = £ a n x n ; e 2d = J2 Pn\ n ; M = £ i n \ n . (3.13) 

n=l n=2 n=4 

The expansion of the 9{ U angles are fixed by Eq.( p.llj) in conjunction with 
the magnitudes of the 9^s in the CKM matrix, as specified in Eq.( [2.2|) p| . 



Having written out the matrices M u , Md according to Eq. (|3.9| ) (with 
Mf», Mf a9 given by Eqs.( |2.9| ) and ( |2.10| )) and expanded each matrix el- 
ement in a "A" expansion similar to Eq.( [3.4|) , one can follow a procedure 
analogous to that in the two-generation case and infer which of the mixing 
angle options give rise to natural mass matrices. The detailed expressions for 
M u and Md are quite lengthy and not tremendously illuminating, so we shall 
omit them and instead only present our findings regarding the naturalness 
conditions. We find that for natural mass patterns, we must require again 
that 

(1) Ou ~ A , 6 lu < A 2 ; 

and, for the other angles, that one of the following options from each set 
holds: 

(2) 9 2u ~ 9 2d ~ A 2 , or 9 2u ~ A 2 » 9 2d , or 6 2d ~ A 2 >> 9 2u ; 

(3) e 3u ~ e 3d ~ a 4 , or e 3u ~ a 4 » e 3d , or e 3d ~ a 4 >> e 3u . 

(3.14) 

The above conditions severely restricts the form of the mass matrices to 
which they apply. As a result, the general expressions for these matrices 
are readily obtained. (The detailed results and their discussion are relegated 
to Appendix A.) Here, as an example, we give a mass pattern with 9 id ~ 
A , 9iu ~ 9 2 d ~ A 2 , 9 2u ~ 9 3u ~ A 4 and 9 3 d ~ A 5 . This pattern corresponds 

5 Notice, in particular, that according to Eqs.(2.4) and ( |2.5|) , S3 = AaX 3 ~ 0(A 4 ), if 
one uses the central values for A and a. 
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to 



M, 



MnA 7 

e iS "u 13 \ 4 

d u X 4 
d 12 \ 3 
e iSd d 13 X 5 



M12A 6 
M22A 4 
M23A 4 

rfi 2 A 3 
rf 22 A 2 
do,X 2 



e~ iSu u 13 X 4 
M23A 4 
1 

e~ i5d d 13 X 5 
d 2 sX 2 
1 



(3.15) 



where the real coefficients Ui/s, di/s are functions of the following 0(1) 
parameters: the CKM parameters A, a; the quark mass ratios £'s; and the 
"A" expansion coefficients {aii, fa, Ts}- The coefficient «i 3 and the phase 5 U 
in addition also depend on the arbitrary phase parameter Sd from the matrix 
A d (Eq. (|310D ) as well as on the CKM phase 5. 

The principal goal of our construction is to allow us to extrapolate LED- 
consistent natural mass patterns to some GUT scale where we can look for 
hints of "new" physics. Nevertheless, low energy (defined here as ~ m t ) 
mass patterns such as the one discussed above are also interesting in their 
own right. For instance, for the pattern given in Eq.( p,15[) , by appropriately 
choosing the signs of the quark masses one can arrange to have 



M, 



(UX 7 + {AafX Q a£ ct A 6 
< ct A 6 £ctA 4 
-Aae iS X 3 



-Aae 




l5 X 3 







V £,db£,sbX 3 
V £,db£,sbX 3 £,sbX 2 AX 2 

AX 2 1 



(3.16) 



with a = {y£,db/£,sb — 1}/A — 0.12. This new pattern now has a large 
number of the much sought-after "texture- zeros" : it has three exact ones to 
begin with, i.e. [MJ23 , [Md]n and [MJ13; and two more to the accuracy 
level 0(A 4 ) of the CKM matrix, i.e. [M u }n , [M u } 12 ~ 0(A 7 ). Moreover, this 
new pattern exhibits useful features commonly exploited in the study of mass 
matrix patterns: sensible CKM and other "predictions" can come about when 
one imposes equalities among matrix elements of approximately the same 
order (e.g. demanding |[M M ]i 3 | = [M u ]22 results in the prediction: \V u b\ — 
m c /m t ), or when one assigns specific values (usually "0") to certain (usually 
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small) matrix elements (e.g. setting [M u ]i2 = results in the prediction: 
sin^c = \Jm d /m s ). 

4 Potentially Successful GUT Scale Mass Patterns 

4-1. A Pattern 

Having constructed certain low energy natural mass patterns, one can 
then apply RGE's to evolve these patterns to some high mass scales where 
global symmetries originating from some GUT texture should become man- 
ifest, hopefully gaining some useful insights. As a study case, in this paper 
we examine the evolution of our natural mass patterns in the MSSM theory. 
For simplicity, we consider only the scenario where the VEV's of the Higgs 
coupled to u-quarks and d-quarks are approximately equal, i.e. tan (3 ~ 0(1). 
The relevant 1-loop RGE's §[|15j are: 

i {(3[/ i[/ ]l 3 -c fc ^)M ii + 3[/i [/ ] i 3M33M3 j } , (4.1) 

l{(3[M33+M4}-4Sfc)M^ 

+ [huh[huh[h D h j } (4.2) 

where the g^s are the three gauge couplings, = (13/15, 3, 16/3), c' k = 
(7/15, 3, 16/3) and hy , h D and h E are the Yukawa coupling matrices for u- 
quarks, d-quarks and leptons,| respectively. To pursue our analysis further, 
we need to solve the above equations to find the mass matrices at the GUT 
scale me — 10 16 (GeV) [@. An input is then necessary at the energy scale m t . 
For definiteness, we assume as a boundary condition [hu{m t )}33 = 1, although 
the general pattern of our result is largely independent of this choice. For the 

6 Although our analysis so far does not concern the lepton mass matrix, one can still 
argue that it must be of the form which reflects the lepton mass hierarchy, especially in 
light of the fact that one wishes to implement the successsful Georgi-Jarlskog GUT mass 
relation p^j : mi,/m T = 3m s /m^ = md/3m e = 1 at some point. It then follows that 
^ r {^i;} ~ [Asl la ~ [^flla > m which case, the contribution of the term tr{h^;} to the 
solution of Eq. ([4.4 ) is very minimal. 



d[h 



U\ij 



dt 

d[h D } i: j 



dt 
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concrete example of our mass pattern of Eq.( |3.15| ), the solution of Eqs. ( |4. 1| ) 
and ( |4.2D gives 



({0.61}un + {0.08}u 2 3 )A 7 {0.61}ui 2 A 6 e-^n^A 4 
M u {m G ) ~ {0.82} | {0.61}ui 2 A 6 {0.61}u 22 A 4 u 23 A 4 

e^ni 3 A 4 u 23 A 4 1 

{0.85}dnA 4 {0.85}d 12 X 3 ({0.85}^^ di 3 + A i3 )A 5 
M d (m G ) ~ {0.38} | {0.85}di 2 A 3 {0.85}d 22 A 2 {0.85}d 23 A 2 + A 23 A 4 

e^di 3 A 5 d 23 A 2 1 

(4.3) 

with A 13 = {0.68}e~ i5u u 13 and A 23 = {0.15}w 23 . The above result illus- 
trates several general features of the RG runnings of natural Hermitian mass 
patterns: 

(1) As is obvious from the form of Md(mc), the Hermiticity of the mass 
matrices is not strictly preserved by the RG evolution. However, the extent 
to which Hermiticity is broken is relatively minor. 

(2) Because of the hierarchy in the mass matrices at m t , the RG runnings 
of various mass matrix elements are quite different. 

(3) This notwithstanding, the likely mass-matrix-element candidates for 
"texture-zeros" at mg are the same ones which are present at m t . These are 
the matrix elements of 0(A 4 ) or smaller. 

These observations suggest a strategy on how to proceed in the search for 
GUT patterns. First of all, since in practice it makes more sense to imag- 
ine the mass matrices at the GUT scale to be Hermitian (or symmetric), 
one should really reverse the procedure. Clearly, if one chooses the GUT 
pattern to be Hermitian (e.g. by manipulating Eq. ( |4.3j ) into its nearest Her- 
mitian form), one should expect only modest deviation in the low energy 
mass matrices from being perfectly LED-consistent. This is even less of a 
real problem since, as a matter of fact, one can obtain LED-consistent mass 
matrices which are non-Hermitian (see the more detailed discussion in Sec. 5). 
Secondly, one can exploit the differences in the RG runnings of mass matrix 
elements to arrange for possible equalities among them at the GUT scale. 
Finally, since "texture-zeros" track between high and low energy scales, one 
can look for possible "texture-zeros" of a GUT pattern in the "texture-zeros" 
or "near-texture-zeros" of its corresponding low energy pattern.[| 

7 Equalities among matrix elements and "texture-zeros' are always desirable in GUT 
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For concreteness, it is useful to demonstrate our ideas with a specific 



example. Let us consider again the mass pattern of Eq.(4.3). By choosing the 
signs of the quark masses appropriately, one finds among other possibilities, 
a potentially "successful" GUT mass pattern in which 



M u {m G ) 



M d {m G ) 




(4.4) 



(4.5) 



where the magnitudes of the parameters are : A ~ 0(1), B ~ 0(A 4 ), C ~ 
0(A 6 ) ; D ~ 0(1), E ~ 0(A 2 ) and F ~ 0(A 3 ). Given this ansatz, one can 
then try to find out whether it fits the LED. As we shall see below, some 
straightforward computation shows that indeed it does.0 



4.2. CKM Predictions 



Running the RGE's backwards, at the energy scale m t the matrices of 
our mass pattern become 



-0.65' 2 


C 


B'e- 1 


a 


1.65' 


B' 


B'e* 


B' 


1 



F' -O.l&B'e-^ 
Si,, ~ I F' 1.2E' 1.2E'-0.18B' 
E' 1 



(4.6) 



Notice that to our accuracy, only Md is slightly non-Hermitian. For ease 
in the computations that follow, it is convenient to specify the approximate 



mass patterns in that they reduce the number of input parameters and, as a result, enhance 
the predictive power of the patterns. 

Alternatively, and in fact more efficiently, we could take as our starting point the 
Hermitian mass pattern in Eq.(3.15) to be our GUT pattern with all the quark mass 
ratios and CKM parameters therein evaluated at the energy scale m,Q (Appendix B) and 
directly arrange for "texture-zeros" and equalities among its matrix elements. 
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magnitudes of the various parameters above by defining B' = bX 4 , C = cA 6 ; 
E' = eX 2 , and F' = /A 3 . With these choices, one can relate the parameters 
b, c, e, / to the quark mass eigenvalues (with the signs of "u" , "c" and "d" 
quark masses chosen to be negative) by solving the corresponding eigenequa- 
tions of the mass matrices, i.e. 

det{M u + £ rf A 4 } = , 
det{M,M]-(^A 2 ) 2 } = ...etc. 

These computations give 

b ~ -0. 66* , 

c ~ ±^(-4.5^ + 0.6^) , 

/ - \J tsbidb ■ 

With these values, one can further calculate the diagonalizing unitary ma- 
trices Ul, Dl from the equations 

ui(M u )u L = M d : a \ 

D{(M d Ml)D L = {M^} 2 ■ 

The final result is 

[CKM] = U\D L 

( 1 Ai 2 A13 \ 

~ -A 12 1 A 23 

V A 31 -A 23 1 / 

where 

Ai 2 = y/U/U A ± + 0-6^ A 2 , 

Ai 3 = 0.7^"^ A 4 ± U^/-4.5U/U + 0.6£ rf A 4 , 

A 23 = £ sb A 2 + 0.7£ rf A 4 , 

A 31 = ^/UUX 3 -0.7^X 4 . 
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Comparing this matrix with the CKM matrix in the Wolfenstein parametriza- 
tion and denoting the absolute values of the quark masses as m g 's, one arrives 
at the following "predictions" : 




raj / m u m r 
sin# c ~ J—±\ - + 0.6— , 

V mc m * 

V cb ~ — + 0.7— , 

m b m t 



V ub * 0.7V±^J-^ + 0^ . (4.7) 
m t m b y m c m 4 

To check the soundness of these results, we choose the "+" sign in the above 
expressions and input various quark mass ratios. Although these are not the 
only possible choices, we find that for 

U ~ 1.55 , U ~ 0.30 ; U ~ 0.73 , U ~ 0.66 , 

we have a decent fit corresponding to the central values of sin^c — 0.22, 
A ~ 0.78 and, ae- i5 ~ 0.31e"^ + 0.05. Choosing = 90° in the last 
equation for example, gives the point (0.05, 0.31) in the p — rj plane, which 
is well within the known constraints (Eq. (|2.6|) ). We note also that the mass 
ratios above are in reasonable agreement with the much more restrictive 
light-quark-mass constraint relation of Eq. (|2.8| ). 



4.3. Other GUT Patterns 

By examining the general expressions for the GUT scale Hermitian ma- 
trices M u and Ma (Appendices A and B), it is not difficult to find other 
potentially interesting mass patterns. We list here four more such patterns 
which are slight variations of the one we discussed in detail above. The first 
two have the same form for the Md matrices as the one in Eq.( |4.5| ), i.e. 



;i&2) M d {m G ) 




9 In terms of the scaling parameter "r" defined in Appendix B, the numerical factors in 
these expressions correspond to r ~ 0.85 (and hence r 2 ~ 0.7, r 3 ~ 0.6). 
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with D ~ 0(1), E ~ 0(A 2 ) and F ~ 0(A 3 ), but have somewhat different 
M„'s: 



1) 



(2) M u (m G ) 








Be -i<P 


f - 


B 


B 


\ Be** 


B 


A 


0(A 7 ); 






1 C 


C 


Be -i<P 


C 


B 


B 




B 


A 



J 



again with A ~ 0(1), B ~ 0(A 4 ) and C ~ 0(A 7 ). The remaining two 
patterns have a different M d matrix which takes the form 

\ 



(3&4) M d (m G ) 
with D ~ 0(1), E ~ 0(A 2 ) and F ~ 0(A 3 ); and the following M u matrices: 



/ F Fe~ 
FEE 
K Fe* ED) 





/ 





\ 


(3) M u (m G ) = 





B 






V o 


B 


^; 


0(A 4 ) and C ~ 0(A 7 ); 








( c 





\ 


(4) M n (m G ) = 








B 




V o 


B 


A) 



with A ~ 0(1), 5 ~ 0(A 2 ) and C ~ 0(A 7 ). The CKM "predictions" of 
the above GUT patterns are most readily obtained in terms of various quark 
mass ratios and the parameter "r" defined in Appendix B, by comparing the 
matrices of these patterns with the general results of Appendix A applied at 
m G (see the example and comments in Appendix B for details). The relevant 
predictions for these GUT patterns are tabulated below in Table 1. 

Several comments are in order at this point: 
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(1) 


(2) 


(3) 


(4) 












sin 9c 


frfid 
V m s 


Jnu m u ^ m c 




fm~d 
V m s 


v cb 


m s 2 m c 

h r 

m b m t 


m s ^ 2 m c 
m b m t 


m s 2 m c 

h r 

m b m t 


V m t m b 


V u b 


m t 


m t m b \ m c m t J 


jm d m s ^_ l4> 

V TO 6 


jm d m s ^_ %<j} 
V m b 



Table 1: GUT pattern predictions for CKM parameters 



(1) Although the mass patterns listed here all contain "texture-zeros", we 
have not tried to impose "texture-zeros" in all possible places. For example, 
one could have in the matrix M u (7tig) of pattern (1) an extra "texture-zero" 
by taking C = 0. The resultant new pattern would, in addition to its CKM 
"predictions", generate a GUT scale quark mass relation corresponding to 
m u m t ~ r 3 m^ which, according to Eq. (|2.9|) , is actually allowed! A systematic 
search for LED-consistent patterns with the maximum number of "texture- 
zeros" has already been thoroughly carried out in Ref. || where, specifically, 
a total of five patterns with five "texture- zeros" were found and discussed in 
substantial detail. 

(2) Because of the specific mass matrix parametrization scheme we have 
chosen, certain frequently-encountered Hermitian mass patterns in the liter- 
ature may not seem transparent from the constructions of our natural mass 
matrices. Still, in general these patterns can be related to our easily derivable 
patterns by some simple unitary transformations. For example, consider the 
following pattern which can easily be arranged from our general results in 
Appendix A, 

/C00\ (OF Ge-^ \ 

M u ~ B ; M d ~ F E E' , 
\ OA/ \ Ge** E' D ) 

where A, D ~ 0(1); E, E' ~ 0(A 2 ); F ~ 0(A 3 ); B, G ~ 0(A 4 ) and C ~ 
0(A 7 ). With the various parameters carefully chosen, this pattern can be 
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transformed into a much more familiar-looking form || 



C 
M' u ~ | O B | : 
A 



by a unitary matrix T (i.e. T M ni d T* ~ M' u ^) with, 

-w 
we** 1 
1 




where u = G/E' ~ 0(A 2 ). 



(3) If one wishes to incorporate the Georgi-Jarlskog mass relation[|TC 



the corresponding lepton mass matrices of the GUT patterns listed in this 
section can be chosen in a straightforward manner. For instance, following 
Ref. B, one can simply let, for patterns 1 & 2, 



(1&2) M,(m G ) 
and similarly for patterns 3 & 4, one can let 



/ F 
F -3E E 
\ ED 





{ ° 


F 


Fe~ l 


(3&4) 


M^ma) = F 


-3E 


E 




\ Fe i(t> 


E 


D 



With these matrices, it is easy to see that the Georgi-Jarlskog relation results 
directly. 



5 Hermiticity Breakdown and Strong CP Complications 

In our introductory discussion we described how to construct Hermitian mass 
matrices from LED information. We must, however, face the fact that imag- 
ining the mass matrices are Hermitian at the weak scale is not a very com- 
pelling assumption. Indeed, as we have argued in the preceding sections it 
is much more sensible to imagine that quark mass matrices are Hermitian 
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(or symmetric) at the GUT scale. When this is the case, the RG evolution 
definitely introduces some non-Hermitian (or non-symmetric) components at 
the weak scale. This was illustrated in explicit detail in the example based 
on the mass matrices given in Eqs.( |4.4| ) and ( f4.5| ). Thus, to be realistic, we 
should instead display a set of natural non-Hermitian weak scale mass ma- 
trices constructed from the LED and then evolve these matrices to the GUT 
scale. 

If one attempts this kind of a general construction from the LED without 
any further constraints, one is immediately faced with considerable arbitrary- 
ness and little progress appears possible. However, if one assumes that the 
resulting weak scale matrices are only "slightly" non-Hermitian, because they 
are Hermitian at the GUT scale, then a general construction becomes feasi- 
ble. In fact, such a construction is really not necessary in the case of SUSY 
GUTs with tan/3 ~ 0(1). In this latter case one can simply detail how the 
CKM parameters and the quark mass ratios evolve to the GUT scale. With 
these parameters in hand one can directly construct natural Hermitian mass 
patterns at the GUT scale. The resulting non-Hermitian mass matrices - by 
construction - will be natural and reproduce the LED. This is basically the 
technique used to deduce Table 1. The details of this procedure is further 
illustrated through an example in Appendix B. 

The presence of non-Hermitian mass matrices at the weak scale, inci- 
dently, raises the issue of strong CP violation. Because of the non-trivial 



nature of the QCD vacuum jl7| , the standard model is augmented by an 



extra CP violating term involving the gluon field strength and its dual 



£> Strong CP = @ F£ V F a ^ v 



The parameter 9 is a linear combination of a phase angle 9 connected with 
the QCD vacuum and another connected with the quark mass matrices 



9 = 9+ Arg{debMi} . 

i=u,d 

One knows, however, that this parameter must be extremely small (9 < 



10~ 9 ) [19], so as to avoid being in conflict with the present bound on the 
neutron electric dipole moment. Why should the QCD vacuum angle be so 
precisely aligned as to cancel (or very nearly cancel) Arg{det M u ^} is not 
known and constitutes the strong CP problem. 
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For the quark mass matrices we have been discussing, if we assume that 
at the GUT scale these matrices are Hermitian then obviously 



Arg{detM u>d (m G )} = . 

However, as we have seen from our analysis, RG evolution induces non- 
Hermiticity. Thus, starting with some Hermitian mass matrices M u ^ at the 
GUT scale, in general, these matrices become slightly non-Hermitian at the 
scale of m t . This is a direct consequence of the RGE's not being Hermitian- 
conjugation invariant. In the 1-loop RGE (Eq. |4.2j ), for example, the term 
[friz] [/%] [/id] is responsible for this non-invariance. Nevertheless, such a term 
is found to be insufficient to generate 

Arg{detM u4 (m t )}^0 . 

In general, however, one expects eventually that at sufficiently high order 
such a term will ensue from the RG evolution. 

The actual order at which a non-zero value for Arg{det M u ^} at the top 
scale appears depends on the underlying theory. For instance, in a globally 
super symmetric theory this mass matrix phase is never generated since it is 
not renormalized |2U|| , while with the standard model it may first appear at 
six loops in the Higgs sector, with an additional gauge boson loop ||21|| . In 
super symmetric theories where SUSY is broken softly the actual contribu- 
tion depends on the breaking. In certain instances no mass matrix phase 
appears |2l| but, in general, if there are non- vanishing elementary or in- 
duced gluino masses one expects a phase to appear For instance, with 
an explicit gluino mass, one induces a non- vanishing Arg{det M u ^{m t )} at 
two- loops P3[ . 

It is quite possible that with the right underlying theory, imposing 
Ar g {det M u ^(mc)} = at the GUT scale suffices to guarantee that 
Ar g {det M u ^(m t )} is much below 1CU 9 . However, this still does not solve 
the strong CP problem unless, somehow, 6 vanishes at mc (which certainly 
is not sufficiently guaranteed by just having Arg{det M u ^(rriG)} = 0). These 
additional observations indicate perhaps compellingly the necessity of having 
some dynamical strong-CP-removal mechanism, conceivably by imposing a 



U(1)pq symmetry |24 
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6 Conclusions 



Interesting patterns of quark masses are surely signals of "new" physics. The 
task of searching for them therefore can be very rewarding. In order to con- 
duct these searches more effectively, we have suggested in this paper the 
idea of natural mass matrices as an organizing principle. This idea, along 
with the efficient mass-matrix-parametrization scheme we have described, 
allows a procedure whereby one can systematically input low energy data 
to construct viable GUT patterns. Encouragingly, this procedure has pro- 
duced a rather small set of "working" mass patterns and our preliminary 
work in extrapolating these patterns to GUT scales has generated some in- 
teresting possiblities. We have discussed, specifically, one such application 
in the context of SUSY GUTs and some potentially successful GUT mass 
patterns were readily found. Although we do not particularly wish to assign 
too much significance to these mass patterns and their predictions, such ex- 
amples do indicate the usefulness of our approach. An important future task 
is to perform a more systematic and complete investigation, with different 
RGE boundary conditions and perhaps different matter contents. 
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A General Expressions for Natural Mass Matrices 

Adopting the parametrization scheme we have developed in Sec. 3.2, we arrive 
at the following general expressions for natural mass matrices, incorporating 
all naturalness requirements (Eq. fl3.14p ): 



M, 



u 



( U\\\ 



unX 6 
U22A 4 
W23A 2 
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Here, 



0(A 9 ) 0(A 8 ) 0(A 6 ) 

0(A 8 ) 0(A 6 ) 0(A 6 ) 

0(A 6 ) 0(A 6 ) 0(A 4 ) 

/ d\\\^ d\ 2 X 3 d\%\^ 



(A.l) 



M, 



d,,\ 3 



doo\ 2 



d,,X 2 



\d* 13 X 4 d 23 X 2 



+ 



/ 0(A 6 ) 
0(A 5 ) 
V 0(A 6 ) 



0(A 5 ) 
0(A 4 ) 
0(A 4 ) 




(A.2) 



Mil 


= U + WlU + \l u e~ i&d - a u A + 


A(e 


«12 


= duict + (3 u [lue~ i5d - a u A + A(e 




Wl3 


= lue~ i5d - a u A + A(e~ iSd - e~ iS ) 


j 


«22 


= tct+ft , 




^23 


Pu 1 




du 


= £,db + C^disb , 




du 


= OLdisb , 




dia 


= lde~ i&d - a d (3dUX , 




d 22 


= £sb , 




d 2 3 


= Pd ■ 





-i8\\2 



}A , 



5 )] 



(A.3) 

The various parameters in the above equations are as follows: X,A,A{= 
aAX -1 ) and 5 are the CKM parameters (Eg. fl2.3| )); £'s are the quark mass 
ratios (Eqs.(|2.9|) and (|2.10|) ); 5d is a free input phase parameter (Eq. fl3.10D ); 
and finally a's, /3's and 7's are input parameters defined from the "A" ex- 
pansion coefficients of the Uj d's (Sec. 3.2): 



sin 9 ld 
sin 9 2d 
sin 9 3d 



a d X 

PdX 2 

ldX A 



sin 6iu = a u X 2 
sin 9 2u = PuX 2 
sin 9 3u = 7„A 4 



The magnitudes of these last parameters are specified below in accordance 
with our naturalness condition (Eq. ( [3.14] )): 



a d ~ 0(1) , a u < 0(1) , with a d - a u X = 1 + 0(X 3 
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< 0(1) , /3 U < 0(1) , with (3 d ~(3 u = A + 0(\ 4 ) ; 
l d < 0(1) , lu < 0(1) , with ld - lu = A + 0(A 8 ) . (A.4) 

The general expressions summarized here are particularly useful for the 
purpose of arranging interesting mass patterns. Specifically, "texture-zeros" 
and equalities among mass matrix elements can, whenever possible, be rather 
conveniently imposed by adjusting the parameters en's, /3's and 7's, subject to 
the constraints in Eq.( |A.4| ). Furthermore, CKM and other "predictions" then 
ensue when the aforementioned constraints overspecify these parameters. 



B Mass Ratios, CKM Parameters and Constructions of Hermi- 
tian Mass Patterns at the GUT Scale 

The relevant formulas for calculating the RG scaling of mass ratios and CKM 
parameters are derived in Ref. [fL5|[25|. Here, we give only a brief summary 
of the results. For the SUSY GUT case we are considering (Sec. 4.1) where 
tan/? ~ 0(1), one finds the following simple RG scaling relations: 

- Mass ratios 

icti.m G ) ~ r 3 £ ci , iut{.m G ) ~ r 3 f ut ; 
isb(m G ) ~r£ sb , idb(m G ) ~ r£ db . (B.l) 

- CKM parameters 

\{m G ) ~ A , 
A[m G ) ~ r A , 

a(m G ) ~ a. (B.2) 
The scaling parameter r in these relations is defined by 

r = e J^I Jo [h uW \ 33 dt ( f = ln ^/ mt }) ( B 3 ) 



which, based on Eq.( |4.1| ) and the 1-loop RGE's for the gauge couplings[j, 
can also be expressed as 



10 



These are, dgi/dt = fo^f /16tt 2 (i = 1,2,3) with b t = (33/5, 1, -3). 
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with 



[hu(m t )} 33 j ~ iV[mGn 



471-2 



[hu{m t )}l 3 I(m G ) 



-1/2 



(B.5) 



V(l*) = Yl{9i(m t )/gi(ii)} 



2ci/bi 



and 



rj(fi) dt 



To obtain a numerical value for r, we input gf{m t )/^ ~ (0.017, 0.033, 0.100) 
as values for the gauge couplings (at writ) Q along with the boundary condition 
[hu{mt)}33 = 1 into the above results, and we find r ~ 0.85.0 

The general expressions for natural Hermitian mass matrices at the GUT 
scale can be gotten by substituting the mass ratios and CKM parameters eval- 
uated at m,Q in Eqs.( |B.l| ) and ( |B.2| ) into the expressions given in Appendix A. 
For illustrative purposes, we "derive" a somewhat generic Hermitian GUT 
pattern and its CKM "predictions" below. 

Choosing in Eqs. (|A.l|) - (|A.3|) the parameters 5d = 5, (3 U ~ 0(A 2 ) and 
demanding u±i, du < 0(X 2 ), one arrives at a mass pattern which can be 
written as 



M, 





C y u Be~^ \ 




f 1 


B x u B . 






x u B A ) 




( 


F VdFe-** \ 




F 


E x d E 


(B.6) 




x d E D j 




0(A 2 ); F 


~ 0(A 3 ); B ~ 0(A 4 ) 


C ~ 0(A 6 ); and 



where A, D ~ 0(1); E 
x's, y's are adjustable parameters which are constrained only by the natu- 
ralness requirement. 



11 These numbers were also used to produce Eq.( |4.3| ). They correspond to a set of values 
for the gauge couplings (at mz) used as inputs in Rcf. Q where, solving the 1-loop RGE's 
with these inputs, th e th ree gauge couplings were found to merge at mc — 10 16 (GeV). 

12 Notice from Eq.( |B.5| ) that the solution for [hu(ma)]33 depends rather sensitively on 
the choice of the boundary condition. However, while still important, this dependence is 
comparatively speaking much milder for r. 
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Mapping the matrix elements in Eq.( |B.6| ) onto those in Eq. ([A.3|) , one can 
immediately establish the following: 

Pu - Xu^ctX 2 , (3 d ~ x d i sb ; 

-a u A + lu e~ iS ~ VuUer iK , 

-adMsb^ + lde~ iS ^ VdadU^e'^ ; 

U + alUX + iyuUfX ~ , 

£d& + al&& ^ . 

Next, solving for the parameters a's, /3's and 7's from the above equations 



and subsequently applying the CKM constraint relations given in Eq. ([A.4|) , 

one has, 



1 ~ v-WU±v-W(^A)-3/2ectA 



±A v /-e itt /(e ct A) - ^ + a*v-6»& A • 

Finally, in the above equations, if one keeps only the significant terms and 
takes into account the RG scaling of mass ratios and CKM parameters (i.e. 
Eqs.( [B.l|) and ( |B.2| )), one sees that the GUT pattern of Eq.( |B.6| ) has the 
following CKM "predictions": 



m d m u m c 

sm6 c ~ J ±J yir 6 — , 

V m s y m c m t 

*cb — -£<f j 



V ub ~ 2/ d — J e 9d -yZr — e ^ ± x d — J y%r 3 — . 

m b y m s m t mj, y m c m t 

The signs of the quark masses above have yet to be chosen as either "+" or 
"— " , depending on which choice is more sensible and gives better agreement 
with experimental measurements of the CKM parameters. 
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